Polchinski equation, reparameterization 
invariance and the derivative expansion 



Jordi Cornelias* 

Departament d'Estructura i Constituents de la Materia 
Facultat de Fisica, Universitat de Barcelona 
Diagonal, 647, 08028 Barcelona, Spain 
comellas@sophia . ecm .ub.es 



Abstract 

The connection between the anomalous dimension and some invariance 
properties of the fixed point actions within exact RG is explored. As an 
application, Polchinski equation at next-to-leading order in the derivative 
expansion is studied. For the Wilson fixed point of the one-component 
scalar theory in three dimensions we obtain the critical exponents rj — 
0.042, v = 0.622 and to = 0.754. 
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1 Introduction 



The derivative expansion Q is, up to now, the most reliable approximation 
available when dealing with exact renormalization group equations (exact RG 
equations, hereafter). Nevertheless, although quite simple at lowest order 
some subtleties appear beyond it. 

They are mostly related to the preservation of the so-called reparameter- 
ization invariance ^| through the derivative expansion. In bold terms the 
problem is the following. It is known that an overall normalization of the ac- 
tion^] does not matter and, therefore, one may define RG transformations in 
such a way that its fixed points (FPs) are normalized to any pre-specified value. 
However, does this freedom survive order by order in the derivative expansion? 
And if not, which is the "best" normalization, in the sense that the expansion 
converges the most rapidly (if at all)? 

The problem is not merely academic because, as we will see, invariance 
of universal quantities under reparameterizations is the ultimate responsible 
of the uniqueness of the anomalous dimension rj, much in the same way as 
the invariance under rcscalings in a linear system over-determines its solution, 
effectively quantizing its eigenvalues. 

These two questions have been successfully answered for equations regarding 
the "quantum effective action" W . Essentially, one chooses RG transformations 
in such a way that some invariance properties of the full equation arc maintained 
order by order in the derivative expansion (see below for details). In contrast, 
for the Polchinski equation J5j, these points seem not even been addressed. 

Specifically, in Ref. || a calculation of critical exponents for one-component 
scalar theories was presented, using Polchinski equation up to second order in 
the derivative expansion. An arbitrary normalization of the action was chosen, 
and a large class of transformations analyzed. Critical indices were found to be, 
not surprisingly, transformation-dependent, as the arbitrary truncation of the 
expansion introduces spurious dependencies. These kind of ambiguities were 
partially solved by the usual criterion of minimal sensitivity J?J , thus choosing 
the transformation which one hopes to converge most rapidly when derivative 
expanded. 

Two problems remain, nevertheless, open. On one hand, as we have already 
mentioned, the possibility of non-standard normalizations was not addressed, 
nor a kind of reparameterization invariance sought. On the other hand, a one- 
parameter family of transformations was left over for which, apparently, a mini- 
mal sensitivity criterion cannot be found. Different transformations within this 
line give slightly different RG eigenvalues although around a FP with signifi- 
cantly different anomalous dimensions. Thus one is left with the impression that 
clear results cannot be obtained and that the derivative expansion as it stands is 
thus of limited use beyond the leading order, at least for those type of equations. 



1 See, for instance, Ref. [jj. 

2 Hamiltonian, in condense matter language. 
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The first problem is specially annoying, because, although we expect that 
different transformations have different convergence properties and, thus, give 
different results, we do not expect that the same RG transformations may give 
a non-unique answer. 

The article is organized as follows. We first review in Section || the deriva- 
tion of Polchinski equation together with its projection. The leading order of 
the derivative expansion is reviewed in Section |^. In Section ^ we explain the 
connection of the anomalous dimension 77 with the reparameterization symme- 
try and its relevance for next-to- leading order computations. A (failed) search 
for RG transformations invariant under a linear representation of the symmetry 
group is discussed in Section |s[ Finally, in Section [| a detailed description 
of our method of calculation, together with the obtained results for critical in- 
dices ry, v (given by the inverse of the relevant RG eigenvalue) and lu (minus the 
first irrelevant RG eigenvalue) of the Ising universality class in three dimensions, 
is reported. Two peripheral subjects are left for the Appendices. 



2 Equation 

Exact RG equations express the change of some quantity (usually the classical 
action) under an infinitesimal RG transformation. Polchinski equation deals 
with one-component bosonic models regulated as 



S[4>] = 77 / P 2 ^~ + SV [</>], (i) 



with some function K(p 2 ) decreasing faster than any polynomial for large p 2 . 
Conventions are as follows. We work on Euclidean momentum space of di- 
mension d, with <pp being the p mode of the Fourier transformed field and 

j p = J-pMp. All variables are dimensionless. Thus in dimensionful units the 
regulating function should be K(p 2 /A 2 ), A = Aoe - * with Ao some fixed scale 
and t the label that parameterize the RG flow. 

RG transformations are usually defined in two steps, first some blocking — 
thinning of degrees of freedom — and afterwards a rescaling of variables so that 
FPs are possible. The first step is accomplished by splitting the field into 



»P-Yp (2) 



with the field 4> p propagated by 

K{p 2 e 2T ) 



(3) 



and <f>p by 



2 



The field </>p^ is then to be integrated out, thus effectively eliminating modes 
mainly within the shell e~ T < \p\ < 1. For infinitesimal r the result is that 
functional integration over the full field (j) p is equivalent to integration over </> p ^ 
with the action 

\ I p 2 K-^e 2 ^^l + S I + r{ f S'^K'ip^- f K'(p 2 )S'l_ p A , 



2 

with 



(5) 

See Ref. |§| for further details. 

Rescaling of variables is carried out with 

p = pe~\ ^=cp p e T ^, (7) 

where r\ is some ^-dependent parameter which coincides with the field anomalous 
dimension in the vicinity of a FP. 
Finally, renaming p — ► p, 

, hm S(t + r)-S(t) 

T >0 T 

= [ S'_ p K'(p 2 )S' p - / K'tf)Sl PiP ~ 2 f p'K-^p^K'ip 2 )^ 

J D J d J r) 



"dp~~ p ' 



dS ^ / cp p S' p - U p p.— S' p , (8) 



where we have (partially) returned back to the full action and written a prime in 
the last term to mean that it serves only to count momenta, thus delta functions 
of momentum conservation should not be derived. 

The derivative expansion amounts to consider a subset of interactions, with 
a maximum power of momenta in their coefficients, and to consider the RG 
Eq. (||) for these terms as if they form a closed subset under the RG, that is, to 
approximate their evolution as if any other operators out of the subset do not 
contribute. 

We deal with the expansion up to two derivatives, with an action 

S[4>] = V[M J d d x + Z[M X -^p 2 K-\p 2 )^ p <t> p , (9) 

which produces two coupled partial differential equations, 

/ = 2K'(0)ff - (fK')f' - (fp 2 K')Z' + d+2 2 ~ 7] f- d ~l + V <pf, 
Z = 2K'(0)fZ' + 4K'{0)f'Z + 2K"(0)f 2 -(fK')Z"-'iK- 1 (0)K'{0)f 
- r] Z- d -^- V Z', (10) 
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with ip = cf>o and f((p) — V'{ip). For convenience, we work with such equations 
after the rescalings 



so that, 



f^ 1 ^ f^hm^ z ^ k ~ 1 ^ en) 



/ = -2ff + f + AZ' + d -±^lf- d -^l^, (12) 
Z = -2fZ' - Af'Z + 2Bf' 2 + Z" + 4/' - r/Z - - - 2 + - 

where the transformation dependencies are bound to the parameters 

, _ (-K>(0))(-fp 2 K>) K»(0)K(p) 
~ (-JK')K(O) ' " (-K'(0W ' 1 j 

These conventions coincide with those of Ref. 0, except for the arbitrariness 
ofX(0).| 

Eigenvalues of the RG control deviations from the FP. They can be cal- 
culated from the linearized version of the RG transformations. That is, by 
taking f((p, t) = f FP (<p) + re xt g(<p), Z(<p, t) = Z FP ((p) + re xt h(tp), with f FP (<p) 
and Z FP (ip) the FP values and r«la coupling constant, 

Xg = -2gf FP - 2f FP g> + g" + Ah' + - + \~ V g - \ + \ g', 
Xh = ABf' FP g' - Ahf' FP - 4Z FP g' - 2h'f FP - 2Z' FP g + h" + ig' (14) 
-rjh - - — y^-<ph\ 
with A being the eigenvalue. 

3 Local potential approximation 

In this section we concisely discuss the results obtained within the first order 
of the derivative expansion — the so-called local potential approximation (LPA, 
hereafter). 

That is, we truncate the action simply to 

s = \JtfK- 1 { P 2 )4>- P 4> P + v{M J d d x. (15) 

The requirement that the non-potential term of Eq. (|l5| ) remains fixed under 



the RG forces r\ = 0, while the evolution of V{f) is (cf. Eq. dl2|)) 



f = -2ff + f + AZ'+ d ^f-^- <pf. (16) 

3 Concretely, after Z -» 1 + 2Z and K(0) = 1 in Eq. <Q), one obtains Eq. (5.1) of Ref. [|. 
The normalization Z(0) = chosen there corresponds thus to Z(0) = 1 in the present paper. 
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The FPs appear as the / = solutions of Eq. (fl6|), 

= -2//' + /" + AZ' + / - ^ tpf. (17) 

Since this is a second order differential equation one may expect, in principle, a 
two-parameter set of FPs. This is not true, of course. First of all, Zi symmetry 
imposes /(0) = 0, thus reducing the arbitrariness to a uniparametric family of 
FPs (labeled by, for instance, 7 = f'(0)). Nevertheless, all but a finite number 
of 7's correspond to solutions which end with a singularity at a finite ip = (po, 

fif) " — as <P fo- ( 18 ) 

Therefore, if we want /(<p) to be defined for the whole range < ip < 00, the 
possible FPs reduce to a finite set. 
The Gaussian FP corresponds to 

f(<P) = 0. (19) 

It is easy to compute its eigenvalues because, after the rescaling tp — 2(pf\Jd — 2, 
the eigenvalue equation, 

g"{p) - 2(pg'{(p) + 2 d+ d 2 J 2 2X g(p) = 0, (20) 

is known to have polynomially bounded solutions, odd in tp (we want to maintain 
Zi symmetry), only for 

X = d-n(d-2), n=l,2,... (21) 

In this case the solutions are the well-known Hermite polynomials of odd de- 
gree 0. The case n = 0, g(p) = corresponds to the identity operator (see 
Appendix |X|). It is curious that this very same eigenvectors appear also within 
the LPA of other RG equations, although the aspect of the exact result may be 
quite different |2j . Being the FP action not a universal quantity, this behavior 
should be probably regarded just as a coincidence. 

One can also identify analytically a FP which must correspond to the non- 
critical high temperature FP, 

f(<p) = V- (22) 
Its eigenvectors are also Hermite polynomials but without relevant directions, 

X = d-n(d + 2), n=l,2,... (23) 

The case n = corresponds again to the identity operator. Because of the sim- 
plicity of the approximation, one can also explicitly write down the renormalized 
trajectory between the Gaussian FP and this one here, 

f(<p,t) = re 2t <p/(l + re 2t ), (24) 
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with r being some (fixed) coupling constant M. 

This latter FP is of no physical interest, but we wanted to quote it mainly 
because it has to be taken into account in the numerical work to follow, in order 
not to confuse it with the Wilson FP. 

Finally, there is also the FP corresponding to the Ising universality class, 
known as Wilson FP. Unfortunately, one cannot obtain it analytically and must 
rely on numerical methods. 

The idea is to integrate Eq. ([l7]) from initial conditions /(0) = and f'(0) = 
7, and scan over different 7's until the correct asymptotic behavior is reached, 

d-2 

f(tp)~ip + Ctp d + 2 +--- as ip — ► 00, (25) 

with C being some arbitrary constant (which one must find out). Note that, 
as explained above, C = coincides with an uninteresting FP which should 
be discarded. For details about the numerical method we refer the interested 
reader to Appendix [b|. 

The eigenvectors are obtained in a similar manner. The linearized RG equa- 
tion at this order is simply 

A.9 = -2gf F p - ZfFpg' + g" + -y-.9 - ^-^W- (26) 

In general its solutions will grow exponentially, 

g(tp) ~ e^ 2 H as ip -> 00, (27) 

whereas for a countable set of A's it will be much smoother, 

d-2-2A 

g(ip) ~ ip d + 2 + ■ ■ ■ as ip — ► 00. (28) 

Choosing the latter behavior, we obtain, from the first two eigenvalues, v = 
0.649 6 and to = 0.655 7 in d = 3. See Refs. §, § . 

Finally, we should mention that with our approach of focusing on f(tp) = 
V'(ip) instead of directly on V(ip) we are always missing the identity operator. 
A discussion about it is postponed to Appendix H. 



4 The anomalous dimension r] 

In this section we explore the connection between the anomalous dimension 77 
and some exact symmetries of the FP action. 

There is a large class of RG transformations which are denoted as linear^] 
because the blocked variables are linearly related to the old ones [|l0|. It in- 
cludes the transformations expressed as exact differential equations. One of 
their main features is that they present always arbitrary parameters, usually in 
their rescaling part. 

4 Not to be confused with linearized transformations around a given FP. 
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These parameters are to be fine-tuned in order to obtain a FP. For instance, 
the Gaussian FP of Polchinski equation, 

^Jp'K- 1 ^-^, (29) 

appears solely after the choice ?y = 0. If 77 is set to some other value (in d = 4), 
either no FP is obtained, or only a non-interacting, infinitely massive FP comes 
out, with no interest for Physics. 

A comment is in order. The operator </>o happens to be the first magnetic 
eigenvector, with eigenvalue Xm = (d + 2 — r/)/2. And, from the known scaling 
relation,! 

(anomalous dimension) = d + 2 — 2Am = 17 = 0. (30) 

Thus the free parameter 77 is directly the anomalous dimension. This is a general 
scenario, although the precise relation between the anomalous dimension and 
the free parameters is not always that simple JlJ, [|. The key point, nonetheless, 
is that there always exists some parameters in the transformations themselves 
which need to be fine-tuned in order to reach a FP, and that these parameters 
are closely related to the field anomalous dimension at the FP. 

On the other hand, any FP is usually associated with a whole family of 
them, with different actions but with the same critical properties. For instance, 
in the Gaussian case above, one finds the line of FPs, 

IJ^K-Hp 2 ) [1 + aK{p 2 )} _1 (31) 

where a is some (real) parameter with aK(p 2 ) > —1 for any p. 

Of course, this is not a surprise, since it is precisely this arbitrariness what 
makes the set of FP solutions finite. It is like the scale invariance of a linear 
eigenvalue problem. The latter may have apparently a solution for each arbitrary 
eigenvalue, but the fact that we can choose the normalization of the eigenvector 
at will over-determines the system, making that only a discrete set of eigenvalues 
are allowed. The fine-tuning procedure, thus, is due to some sort of underlying 
arbitrariness (some sort of symmetry). Had a FP equation no invariances at all, 
no parameter would have to be fine-tuned. 

The form of the eigenperturbations about the FP also vary along the line of 
FP. For instance, the mass term is 



+ aK{p 2 )}- 2 ^ p cj> p , (32) 
with eigenvalue 2. 

Therefore, to summarize, linear RG transformations present two important 
features: They contain free parameters which have to be fine-tuned in order a 
FP to be reached; this feature reflects an underlying symmetry which manifest 
itself in a line of equivalent FPs. 

5 See, for instance, Rcf. 
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So far for exact computations. Let us now discuss the actual case, that is, let 
us discuss how this picture survives under the (up to now) unavoidable approx- 
imations of RG equations. For definiteness, we take the derivative expansion 
at second order, Eq. (p"2|), but the following should be easily extendible to any 
feasible truncation. 

First of all, it is clear that the line of FPs, Eq. ( |3l| ) , is hard to reproduce well: 
Each FP has different 0(p 4 ) terms and higher and only those with such terms 
comparatively small will be correctly described by our truncation. Therefore, 
one immediate consequence is that it is not that harmless to blindly choose 
one normalization (one value of a) of our FP action when some truncation is 
involved; e.g., in the Gaussian FP above, several choices of the parameter a will 
be quite faithful while others will probably poorly capture the properties of the 
FP. 

Moreover, the opposite phenomenon also takes place: One also finds trun- 
cated FPs for the wrong choice of the transformation parameter 77. To un- 
derstand better this annoying feature, let us recall Wegner-Houghton equa- 
tion||§, 




A line of FPs is obtained here with 

J\ P \<i 

and 77 = e. This phenomenon is of quite different nature as the similar one 
described above. There the whole line shares the same critical properties, here 
it does not; there we have well-behaved actions throughout the fixed line, here 
nearly all of them are terribly non-local (in the sense that we cannot expand 
the action integrand in a power series of p 2 ). What happens is that we have one 
physical FP (the e = case) and a line of spurious ones. 

For exact calculations this is not a problem since RG transformations al- 
ways map local actions into local actions and thus only the e = FP may be 
reached. Nevertheless, simply studying V(ip) and Z((p), one may not be able to 
distinguish non-local FPs from the local, physically meaningful, one. 

The recipe to cope these problems is not evident. We have, nevertheless, 
some hints. Clearly, one should succeed if one manages to compute a FP as local 
as possible (ideally, one FP with only the V(ip) and Z(ip) terms). The problem 
is to try to reach this goal without computing the next order of ones expansion. 
The key is to try to check, within the given order, some known property of the 
exact solution, and take the approximate FP which best reproduces it. 

A glance at linear systems again may help. Let us imagine that in an eigen- 
value problem of linear algebra, one is so bold to approximate it in such a way 
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that eigenvalues are no longer invariant under re-normalizations of the corre- 
sponding eigenvectors. One would probably find no quantization of eigenvalues. 
However, if one insists, after all, to make sense of the approximation, a clever 
choice would be to look for the eigenvalues which are, to some extent, invariant 
under rescalings of their eigenvectors. This is, in fact, the method used from 
time to time for exact RG equations, following the procedure of Ref. || 0]. 

The trick is to find the FP which is first-order (at least) independent of the 
free parameters. This feature results in the presence of a redundant operator Q 
which is marginal or nearly marginal. It has to be redundant because it may not 
generate true RG trajectories, but simply reparameterizations of the theory; it 
has to be nearly marginal because it generates, at least approximately, a line of 
equivalent FPs (the equivalence being guaranteed by the redundancy). 

In our computations, therefore, we first plot the anomalous dimension versus 
the normalization of the kinetic term. It is seen that this curve presents always a 
maximum, which we regard as the best approximation for the actual anomalous 
dimension. This is further checked by the finding of a nearly marginal operator, 
which is known not to exist in the physical spectrum. See Section ^ for the 
details. 

Before ending, let us make one further remark. The criterion explained so 
far relies on an invariance property of the physical, local, action. This does not 
mean that non-local FPs do not present their own set of symmetries under repa- 
rameterizations. But the claim is that the properties of local solutions should be 
well-described by local approximations. On the contrary, the analogous of the 
redundant marginal operator found for local FPs should be terribly non-local 
and, thus, poorly mapped within the truncation. So poorly that they do not 
even show up S. 



5 Linearly symmetric RG transformations 

Before going on to actual computations, we would like to address one more issue. 
One knows that the accurate anomalous dimension r\ should be close to the one 
for which the action presents an (approximate) reparameterization symmetry. 
That is, why do we not choose RG transformations for which the symmetry is 
linearly realized, and try to keep this realization through the derivative expan- 
sion? This is, in fact, the approach taken in Ref. Q for the effective-action type 
of exact RG equations. 

The search for such class of transformations is reported in this section. The 
conclusions turn out to be, unfortunately, rather deceptive. We find a class of 
regulators with a linearly realized reparameterization symmetry for the exact 
equations; but, first, the symmetry is broken at finite order in the derivative 
expansion; and, second, the regulators associated with such transformations 
turn out not to regulate, at least not in a finite order in the derivative expansion.^ 
Therefore, one is lead to the general scenario outlined in the previous section: 

6 See |15[. and references therein. 
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One must scan a range of normalizations and decide on some kind of locality 
criteria which is the most reliable one. 

Let us seek some functions JC(p 2 ) and T(p 2 ) such that the transformations 

p = Xp, K(\ 2 p 2 )=JC(p 2 ), fa, = (35) 

on an action converts it into another action that satisfies the same RG equation.^ 
In order to obtain an action that satisfies again a Polchinski-kind equation, 
these functions must satisfy 

Kr x K! -T- X T' = K~ X K' 

\- d - 2 T- 2 K! = K', (36) 

with solutions 

for some function a(A) and a real number n. 

If we further require the RG transformations to be exactly the same as those 
of the initial action, the regulator must satisfy 

K{X P } - l + a(X)K(p 2 y (38) 



K (P 2 ) = ^^72 ' (39) 



which implies 

n/a 

i + {fjey 

with a = — a'(l) and q 2 , n real non-negative numbers. This are the regulators 
previously reported in Ref. 

If we were to use the above regulators, then, for every FP, there would exist 
a whole line of them, corresponding to different parameterizations of the action 
(different choices of the dummy variables <f> and p). The FPs within the line 
would be related by simple linear transformations, generated by a marginal re- 
dundant operator, as corresponds to any reparameterization. All the FPs would 
be physically equivalent as long as infra-red properties are concerned jli], ||. 
And if these RG transformations were used, then the discussion on normaliza- 
tions would be totally void. Unfortunately, as we have pointed out above, there 
are important objections. 

First, this kind of symmetry does not in general hold order by order in the 
derivative expansion, unless a(A) = 0, as can be easily stated by considering 
Eq. ©. 

Second, the above class of functions K{p 2 ) do not regulate, at least they 
do not when we consider a fixed order in the derivative expansion, because 

7 One may try more general scenarios, like T = F(p, 4>p)- Nevertheless, not even in these 
general frameworks, any further invariant transformations have been found besides the ones 
considered in the text. 
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Figure 1: Feynman diagram representing the second term in Polchinski equation. 
The shadow area stands for a generic vertex (which contributes with a coefficient 
of 0(p 2n ) at n-th order in the derivative expansion and the thick solid line stands 
for a propagator oc K'(p 2 ). As usual one has to integrate over the closed loop. 

they behave asymptotically as polynomials. That is, a simple study of one-loop 
diagrams (see Fig. fil) show, for instance, that the second term of the RG Eq. (p) 
is not finite at 0(jr") in the derivative expansion with regulators of Eq. ([39) 
whenever n > n + 1 — d/2. 

The conclusion is, therefore, that the symmetry cannot be linearly imple- 
mented order by order in the derivative expansion, in contrast to other equa- 
tions [Q for which the maintenance of reparameterization invariance selects a 
preferred scheme. Polchinski equation in this sense is unfortunate: One cannot 
avoid scanning over different normalizations and selecting a criterion to dis- 
criminate among them. The true FP corresponds only to one of the possible 
normalizations; the rest should probably describe non-local actions. 

Finally, let us look again at the Gaussian FP, Eq. (|3l|). It can be found 
with a regulator of the form of Eq. ([?9|), and in such a case we have a nice 
example of the symmetry we have talked about. Our second objection does not 
apply here, since the Gaussian FP is ultra-violet finite anyhow. Nonetheless, 
the example shows that linear transformations are not the whole story, since, 
even for arbitrary choices of the regulating function K(p 2 ), the line of FPs exist 
and it has to be generated by non-linear transformations. This last possibility 
is what we are seeking for the Wilson FP in d — 3. 

6 Results 

6.1 Gaussian FP 

We have already solved the FP equation within the Gaussian FP exactly, and 
also for the LPA. Let us briefly review how it appears in second order in the 
derivative expansion. 

It consists solely on a kinetic term with a constant potential, that is, 
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And, from the second equation in (]12j), r\ = 0. 

The eigenvectors have a quite simple expression in the case a = 0. Namely, 
there appear the same Hermite polynomials for g((p), with h(tp) = 0, obtained 
when dealing with the LPA, with eigenvalues A = d — n(d — 2). 

There are also other eigenvectors with (for a = 0) g(ip) = and the func- 
tion h(ip) a solution of 

2A 

h'\(p)-2yh\0)-2-^h{$)=Q, (41) 

with (p = v/ '^ 2 (p. Again they are Hermite polynomials, this time of even degree, 
and their eigenvalues 

\ = -n(d-2), n = 0, 1,... (42) 

The general form for these last perturbations of the FP action is, for the 
n-th such operator, 

i /d*(d0) 2 2n (z) + ---, (43) 



with the dots standing for terms with less number of fields. All of them are 
irrelevant but the one with n — 0, which is a marginal redundant operator, 
\ J dx{d(j)) 2 , that generates a line of equivalent FPs. It clearly corresponds to 
a change of normalization of the original FP action. 

The general case a ^ is technically a bit more complicated. The eigenvalue 
equations are now 

Xg = g» + Ah' + ^g-^<pg', (44) 

aK(0) , ,„ d-2 
l + aif(0) y 2 ^ 

The marginal redundant operator is still g(ip) = 0, h(ip) = constant, of course. 
And the mass operator also conserves its form, g(cp) = ip, h{ip) = 0, but it is 
clear that the tower of eigenvectors of the form of Hermite polynomials will not 
survive and a more complicate computation has to be done. For instance, the 
eigenvector with eigenvalue A = 4 — d (the "J d d x 4 (x)" operator) is, for d ^ 4, 

m = «f-("+"i^5y)*. <«> 



/i(^) = 2 



aif(0) 
1 + aX(0) 



6^-3(1+ ^ ^ (0) 



d 1 + aif(0) 



Clearly, the FP with a = plays a special role in the derivative expansion as 
far as convergence properties are concerned. 

We do not discuss it any longer here, as we want only to look at it as an 
illustration of the kind of problems that might appear for the, more interesting, 
Wilson FP. 
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6.2 High temperature FP 



The analysis of the high temperature FP is also easily extendible to one further 
order in the derivative expansion. It appears as f(ip) — ip, as before, and 
Z(ip) = Z, with Z a constant related with rj by the equation 

C2 - n) 2 - - 

0=^ — Y~ B + 2 ( 2 - i) - 2 ( 2 - v) z - VZ- (46) 

Note that now r\ is totally arbitrary and there is no evidence of a line of 
equivalent FPs.^j Furthermore, for some choice of r\ we have Z = 0, thus making 
self-evident its physical interpretation of an infinitely massive theory with zero 
correlation length.^] If Z ^ 0, the same physical interpretation shows up when 
considering dimensionful variables, 



a 2 r , z 

~2 



r"P + -77 I P 2 <P-p<f>p H , (47) 



and letting A — > oo. 



6.3 Ising universality class: FP 

We now discuss the Wilson FP in d = 3. First we sketch the numerical method to 
find it, leaving the details for Appendix |b| Then we turn to the normalization 
suitable for dealing with the reminiscence of the line of FPs present in the 
exact computation. Finally, we asked ourselves about the best choice of RG 
transformations. The answer is given in terms of rapidity of convergence of the 
derivative expansion, proposing a criterium which extends that of Ref . || . 
The FP action is the solution of RG Eqs. (||) with f(ip, t) = Z(cp, t) = 0, 

= -2ff + f» + AZ> + d -±^f- d -^Lpf, 
= -2fZ' -AfZ + 2Bf 2 + Z" + Af 

-nZ-^-^-cpZ', (48) 

and the boundary conditions 

/(0) = Z'(0) = 0, 

f(<p) ~ + Cip d + 2 -i H , as ip — > oo, (49) 

Z{tp) ~ D -\ , as ip — > oo, 

where C and D are constants (which have to be determined). The first two 
conditions come from imposing ^-symmetry, while the last two come directly 



from the FP Eqs. (48), once we require the solutions to exist for the whole range 
< (p < oo. 



8 We mean FPs with the same, e.g., rj. 
9 See, for instance, Ref. IfLql. 



13 



That is, from the asymptotic conditions we have three free parameters (C, D, 
77). The use of Z2 symmetry imposes two constraints. Only after fixing one more 
condition, e.g., Z(0) = 1, the system is completely determined. One obtains, 
thus, (apart from the two special FPs discussed above) one approximation for 
the Wilson FP for any normalization Z(0). This is to be compared with the case 
where one manages to maintain explicitly a reparameterization symmetry 
In the latter case the invariance allows to arbitrary fix the remaining parameter 
to whatever value we choose, just as in a linear system one is allowed to fix 
the scale of the solutions. On the contrary, our results depend on this last 
parameter Z(0). 

In Figs. H [3] we plot two examples. The first one is for A = 0.53, B = 0.40, 
which are the values that, as explain below, we take as the most reliable. The 
second plot is for the values A — 0.8, B — 0.5 which are the ones taken in 
Ref. 0. Note that the normalization chosen in that reference (Z(0) = 1, in 
our notation) has nothing special from our point of view, while the appropriate 
normalization for which reparameterization invariance is (partially) recovered is 
Z(0) = 1.25, with 7] = 0.51. 

As explained earlier, this pattern is a consequence of the breaking of repa- 
rameterization symmetry. The line of equivalent FPs is converted into a line 
of non-equivalent FPs because only the most local ones can be adequately de- 
scribed by our approximation, being the rest quite distorted. Furthermore, 
the latter can be easily confused with some approximation of genuine non-local 
ones. And, as explained above, we take as a hint for locality the fact that a 
reminiscence of the line of FPs is present. That is, we take as the closest to 
the exact Wilson FP the one with an anomalous dimension which is first-order 
independent of Z(Q). For a given transformation (for fixed A and B) there is 
only one point with that property. 

Once we know how to compute the anomalous dimension given the RG trans- 
formations, we have to ask ourselves which are the most appropriate blockings. 
This is an old subject, which has been studied mainly on the lattice. For ex- 
act RG equations and their approximations it was addressed in the pioneering 
work of Rcf . H . After it no much progress was achieved until Ref. [|| , where 
a systematic search over different transformations was made and a proposal of 
"best" transformation was given, based on a reformulation of the principle of 
minimal sensitivity used in perturbation theory. 

Unfortunately, this latter reference is incomplete in two important direc- 
tions, as we have already outlined in the introduction. The first one is its poor 
calculation of the anomalous dimension, without taking into account the prob- 
lem of the breaking of reparameterization invariance of the theory, thus making 
their results not quite trustworthy. The second problem is that their analysis 
did only half the way, leaving rj strongly dependent of one free parameter. 

Nevertheless, we think the analysis of Ref. || is basically correct, as far 
as the general dependencies on A and B are concerned. That is, if a three 
dimensional plot of 77 as a function of A and B is made, one identifies a sort 
of "hollow" around which rj presents a one-dimensional minimum whenever the 
hollow is crossed. To plot it, it is best to fixed one parameter, say B, and plot 
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Figure 2: Plot of the anomalous dimension rj for various FPs, with fixed 
A = 0.53, B = 0.40, as a function of their normalization Z(0). These are 
the transformations we use later on to calculate eigenvalues. The maximum of 
the curve occurs at Z(0) ~ 1.31 with r\ = 0.042. 

r\ as a function of the other one, A in our case. These kind of plots are shown 
in Fig. ||, whereas a sketch of the location of the hollow in the A-B plane is 
presented in Fig. ||. 

We choose, as in Ref. [||, a minimal sensitivity criterion to select the trans- 
formations for which r\ yields on the hollow. This should be the ones for which 
the derivative expansion converges the fastest. See Ref. |?J for an introduction 
of minimal sensitivity ideas. 

Nevertheless, the strongest dependences of rj on the transformations occur 
along the hollow, not across it. Compare, for instance, a plot of different 77's 
along the hollow (Fig. ^) with Fig. [| above. In fact, the reason which prevents 
us to 3d-plot r\ versus the A-B plane is that the quite different scales make it 
difficult to visually appreciate the details of the hollow. 

We must then establish a criterion in order to keep only a small range of 
B values if we wish to have a reasonable prediction for 77. If not, it can take 
nearly any value we want, since, as it was already observed in Ref. ||, it grows 
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Figure 3: Plot of the anomalous dimension r\ for various FPs, with fixed A = 0.8, 
B = 0.5, as a function of their normalization Z(0). These are the transforma- 
tions used in Ref. ||. The maximum of the curve occurs at Z(0) = 1.25 with 
T) = 0.051. 

nearly linearly along the hollow. The other parameter one has to fine-tune to 
obtain the FP, e.g., 7 = /'(0), behaves also in a similar fashion (see Fig. fj]). 
For consistency with previous choices, the criterion we choose must be based in 
convergence properties of the derivative expansion. Clearly, one cannot focus on 
77, since it changes from 77 = 0, imposed by the LPA, to a value close (hopefully) 
to the actual one. This is because it has to do with the newly added term, 
Z((p), which was not present previously. However, if the expansion converges 
sufficiently fast, the function f(tp) obtained at second order must not be too 
different from the one obtained at lowest order. This is the criterion we propose: 
One should trust those FPs whose values of 7 are close to the same, unique, 
number obtained from the LPA. This selects the transformations around those 
with A = 0.53 and B = 0.40, with 7 = -0.229 (to be compared with 7 = 
—0.228 6 from the LPA). Its anomalous dimension is 77 = 0.042. The FP action 
is plotted in Fig. ||. 

Note that with this choice we hope to reproduce well the behavior of the 
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Figure 4: The anomalous dimension rj as a function of A for B — 0.30 (solid 
line), B = 0.40 (dashed line) and B = 0.50 (dotted line). 

FP action near <p = 0. We do so following the common lore that this is the 
relevant part to reproduce if quantum fluctuations of the field are not very 
important, i.e., whenever rj is small. In fact, the whole idea would fail if the 
LPA was not sensible enough, which in turn is based on t] = not being so 
crude an approximation. It is noteworthy, nonetheless, that there are other 
transformations which reproduce better than ours the asymptotic properties of 
f(<p) but give quite worse results. For instance, in Fig. ^ we plot the function 
f(<p) for the LPA, together with the chosen A = 0.53, B = 0.40 FP and with 
that from A = 0.21, B = 0.70. The last FP coincides quite well with the LPA 
one for ip — > oo but differs from it significantly for ip ~ 0, just the opposite 
behavior of ours. Nevertheless, our choice seem to give a better anomalous 
dimension (77 = 0.042 instead of i] = 0.058). 

Before finishing this section, we would like to insist that we have just pre- 
sented a standard calculation within the derivative expansion up to second order, 
without any further approximation. The discussion above is focused on the set 
of transformations best suited for such a calculation. Without a clever choice, 
the expansion is poorly convergent and higher orders are necessary to provide 
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Figure 5: The "hollow" of rj in the A-B plane. 



the desired accuracy. 



6.4 Ising universality class: critical indices 

The computation of the eigenvalues is similar to that of the anomalous dimension 
•q. We have two asymptotic free parameters g\ and h\ for a polynomially growing 
g(<p) and h(<p), 



h(<p) 



9 if) ~ 9i<P d+2 ~ r > + ■ ■ ■ as ip — ► oo, 

-4+2tj-2A -8 + 277-2A 

fiQ(p d + 2 ~v + h\ip d + 2 -n + ■ ■ ■ as ip — > oo, 



(50) 



plus the eigenvalue A (ho is a parameter directly determined from gi), but 
Z 2 symmetry imposes g(0) = h'(Q) = and the linearity of the eigenvalue 
Eq. ( prf ) allows us to fix on more parameter, say g'(0) = 1. In this way the 
system gets determined and we encounter only a countable number of solutions. 
We discuss the three most relevant ones.[^] We concentrate here on the FP with 
A = 0.53, B = 0.40, d — 3, as we stated above. 

10 Recall that there is also the hidden solution A = 3. See Appendix 
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Figure 6: The plot of 77 along the hollow, parameterized by its B value. This 
plot is to be compared with that of Fig. 0. Here the curve is much steeper. 



The only relevant eigenvalue found gives v = 0.622. 

The second eigenvalue is not part of the physical spectrum, but corresponds 
to the redundant operator which generates the line of FPs in the exact case. In 
fact, its appearance is a crucial test of consistency of our method for calculating 
rj. This is so because it must necessarily be marginal, otherwise it cannot connect 
two points which both of them are FPs of the same RG transformations. For 
exactly A = .53, B = 0.40 and Z(0) = 1.31 we obtain A = 0.001. A plot of this 
redundant eigenvalue for different ?7's is given in Fig. [l0|. Note that for each 77 
there are two values of A, corresponding to the two FPs which exist (cf. Figs. |[ 
||). The eigenvalue vanishes only when there is a unique FP given 77, i.e., when 
77 is approximately invariant under Z(fi). 

Finally, we seek for the first irrelevant eigenvalue, which controls the first 
deviations to scaling.^] It gives ui = 0.754. 

Together with the eigenvalues, there comes also their corresponding eigen- 
vectors. As an example, we plot in Fig. [lljthe one corresponding to the marginal 
redundant operator. 

11 See again, for instance, 
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Figure 7: The plot of 7 along the hollow, parameterized by its B value. 
6.5 Accuracy 

All figures given so far are significant. Such a precision is achieved because, 
ones the problem has been set up, one can push the numerical calculations 
sufficiently far to obtain a high accuracy (of even five or six digits). In fact, the 
main sources of errors are, on one hand, the problem introduced by different 
convergence properties of different transformations; and on the other, the poor 
knowledge of the corrections one may expect at higher orders. 

The first problem can be turn to our favor. Once our criteria is accepted, 
there is only left the problem of to which accuracy we want 7 from the LPA 
be the same as the calculated at second order. Clearly, we expect the quantity 
to be corrected somehow at each order in the derivative expansion, so it is 
not sensible to ask for a terribly high coincidence. If we restrict ourselves to 
FPs with 7 three significant digits equal to the LPA value, then the anomalous 
dimension varies within the range -q = 0.041 9-0.042 6, thus providing us with 
an indication of the accuracy of our results. The eigenvalues arc, nonetheless, 
far more insensitive, specially the relevant one. As an extreme example take, 
for instance, A = 0.21, B = 0.70, which gives v = 0.631, u> = 0.689. 
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Figure 8: The FP functions f(ip) (solid line) and Z(ip) (dashed line) for second 
order in the derivative expansion with A = 0.53, B = 0.40 and Z(0) — 1.31. 

The main drawback of the whole framework, nevertheless, is its poor treat- 
ment of systematic errors introduced by the truncation of the expansion. In 
fact, there are, up to now, no methods to estimate the contribution from higher 
orders. It is amazing that the same feature that makes the technique so pow- 
erful (its no use of a small parameter like e = 4 — d or the inverse number 
of fields 1/iV), also prevents us to make a thorough study of its errors, thus 
reducing the reliability of its conclusions. 

Exact RG in this respect is special: in contradiction with what usually occurs, 
it happens to have the peculiarity that it is computationally quite simple, but 
its restrictions arrive in its deficient treatment of errors. We have not a clear 
way to control the size of next corrections within the derivative expansion, nor 
wee can evaluate how good our criterium for improving convergence is. This is 
the main problem of the method, not the computation itself.Fj 



12 Gauge theories are an exception. In this case the main problem is to find a non- 
perturbative gauge-invariant expansion to deal with the RG equation [11 71] - 
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Figure 9: Comparison of the function /((/?) from the LPA (solid line) with the 
second order ones corresponding to (1) A = 0.21, B = 0.70 (dashed line), which 
becomes similar to the previous one for large tp, but differs significantly from it 
for ip ~ 0; (2) A = 0.53, B = 0.40 (dotted line), which do not match with the 
first one for large tp, but coincides with it quite well for ip ~ 0. 



6.6 Discussion 

Table | contains a summary of our results, together with a comparison with 
results from the LPA and also from other methods (exact RG for the effective 
action and a combination of best known estimates M). It is remarkable that 
with exact RG computations one may obtain numbers quite close to the best 
known estimates, and with considerably less effort. 

We hope that the present paper serves to convince the reader that reported 
ambiguities in universal quantities computed with Polchinski equation || [H| 
do not constitute a problem, but rather can be turned to our favor .p| There is 
still open, nonetheless, the problem of how one can compute systematic errors, 
which is the reason that prevents exact RG methods to be totally competitive 
with the best existing techniques (Monte Carlo RG jl9| , e-expansion [^0| , . . . ) . 

13 See also Ref. Jl5). 
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Figure 10: The redundant eigenvalue A as a function of 77 for A = 0.53, B = 0.40 
and varying Z(0). 

As a final comment we want to recall that exact RG methods seem to give 
always an anomalous dimension 77 that increases from zero (LPA) to a value 
which is slightly above the best estimate; while the critical exponent v goes the 
other way round, going form a high value from the LPA to another one which is 
probably a bit low. The exponent w does not seem to oscillate around its actual 
value, but this may only reflect our poor knowledge of it. The apparent general 
behavior of r\ and v is, however, quite intriguing, because if it is really general 
then it may constitute a first step towards an accurate estimation of errors. 
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Figure 11: The marginal redundant eigenoperator for A — 0.53, B — 0.40 and 
Z(0) = 1.31. The solid line corresponds to the function /((/?), while the dashed 
one to Z(ip). For display purposes, we have chosen the normalization h(0) = 1, 
which gives g'(0) = -0.031. 

A Identity operator 

It is known that, at any FP, there always exists an eigenoperator with eigenvalue 
A = d, which is the highest dimensional operator of the theory.^ Let us briefly 
review the arguments that support it |plf . 

Let us take any FP action S* and consider perturbations with a set of 
operators O a , 

S = S*+j a O a , (51) 

where j a are some coupling constants. Implicit summation over repeated labels 
is understood throughout. 

14 Or the lowest dimensional operator in particle-physics language. Recall that in condensc- 
matter language one usually talks about the dimension Dqm °f an operator to refer to the 
(quantum) dimension of its coupling (in units of mass). This is natural because these are the 
fields which appear in the partition function of the theory after the path-integration. Particle 
physicists, however, usually talk about the dimension Dpp of the operator density itself, 
which is, trivially, Dpp = d — Dqm ■ The identity operator thus have Dpp = 0. 
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LPA 


Polchinski 


eff. action 


best known 


V 





0.042 


0.054 


0.035(3) 


V 


0.650 


0.622 


0.618 


0.631(2) 


UJ 


0.656 


0.754 


0.897 


0.80(4) 



Table 1: The critical exponents 77, v and w for (1) the LPA of Polchinski equa- 
tion; (2) derivative expansion at second order of Polchinski equation; (3) deriva- 
tive expansion at second order of the effective action RG equation [EJ ; (4) com- 
bination of best known estimates taken from Ref . W . 



The partition function is then 

Z= f 2?0e- s *-^°=, (52) 



and we define the thermodynamic densities 

Ma = jj ttt- In Z , (53) 

with V the volume of the system (needed in order M a to be an intensive quantity 
and, thus, defined in the thermodynamic limit). 

After a RG transformation, the system is shrunk in all linear dimensions and 
new couplings j'p are introduced, 

V' = e- dT V 7 #=j£0«), (54) 

with the partition function kept invariant. The above defined densities thus 
transform to 

M a = e-*"^ A InZ = , "'" A.i/' (55 ) 
Around a FP, the linearized RG transformations are precisely the matrix 

di 1 

■x = ~dj~~ ' w ith the 
at a FP, M' = Mp, 



lZ/3 a = ^f-, with the right hand side evaluated at the FP. Therefore, using that, 



n f3a M fj = e dT M a . (56) 

The final result is, therefore, that the densities {M a } form an eigenvector of the 
linearized RG transformations, with eigenvalue A = <i.p| 

However, the above densities are expectation values of local operators and 
we may worry if all of them vanish. That this is not the case is due to the 
existence of the identity operator, whose expectation value is 

<l>=iAln [ T>d>e- s ' +joV = 1. (57) 
Vdjo J 



15 Again, the terminology is somewhat confusing. We have usually talked about eigenvectors 
to mean eigenoperators, that is, O = M a O a in this case. 
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This completes the proof. 

With Polchinski equation, one may also trivially show that the identity op- 
erator is always an eigenvector with eigenvalue A = d (just substitute S = 
S* + J d d x in Eq. (§) and use that S* is a FP). 

In the derivative expansion we do not notice it because it is g{ip) = h(ip) = 0. 
However, if the FP potential is Vfp{<p) and its eigenperturbations v((p), then 
the linearized RG equation for the LPA may be rewritten as (cf. Eq. |2^) 

Xv = -2v'V FP + v" + dv- ^T^W, (58) 
with X = d, v = constant, a trivial solution. At next order we have to solve 

Xv = -2v'V FP + v" + Ah + dv - d - 2 + 71 tpv', 
Xh = -2V FP ti - 2v'Z' FP ~ AVpph - Av"Z FP + <iBV FP v" (59) 
+ h" + 4v" - V h- d - 2 ^ 

And, again, a trivial solution is 

A = d } v = constant, h = 0. (60) 



B Numerical methods 

We explain in this appendix the numerical methods used at second order. The 
calculation of the LPA results can be inferred form them and, moreover, they 
are discussed elsewhere ||. 

The FP is searched by shooting from the origin to some point ipo where we 
impose the asymptotic conditions valid for tp — > oo, (cf. Eq. (ff9|)), 

f M - V d-2 + V 

fO— " 2^ = -T—7, -> Z (W = 0. (61) 

Of course, these conditions cannot be imposed for finite ipo, but we turn this 
fact to our advantage. We move ipo between 4 and 6, calculating r\ and 7 at 
each step and we observed that the results stabilizes for the five or six significant 
figures at some point between 5 and 6. Thus it serves not only to check that 
the sub-leading asymptotic terms can be effectively discarded but also to check 
the accuracy of the overall numerical method. 

Incidentally, all numbers quoted are obtained always with at least one signif- 
icant digit more than shown, to make sure that the numerical errors are under 
control. 

The computation is made with one normalization at a time, increasing it 
with finite amounts of 0.01. From the FPs so obtained, we keep the one with 
maximum 77, as explained in Section |^. We check that this grid does not affect 
the results, within our accuracy. For instance, for the quoted ij = 0.042, A — 
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A 


B 


Z(0) 


f(0) 


V 


0.55 


0.39 


1.31 


-0.23017 


0.0413 


0.53 


0.40 


1.31 


-0.228 99 


0.042 


0.52 


0.41 


1.32 


-0.228 07 


0.042 6 


0.49 


0.42 


1.33 


-0.22711 


0.043 3 



Table 2: Comparison of different FPs on the "hollow" (see Section ||). Recall 
that /'(0) = -0.288 60 within the LPA. 

0.53, B = 0.40, we have: for Z(0) = 0.30, rj = 0.041966; for Z(0) = 0.31, 
r) = 0.041 985; and for Z(0) = 0.32, 77 = 0.041 981. 

The scanning of the A-B plane is made in a similar fashion (also with finite 
amounts of 0.01). To identify the hollow mentioned in Section || we fix first B 
and obtain the FPs for different A, and choose the one with the minimum 77 (see 
Fig. H). Again, the grid should not be a great concern. Compare, for B = 0.40: 
A = 0.52, 77 = 0.041 985 1; A = 0.53, 77 = 0.041 984 8; A = 0.54, 77 = 0.041 989 1. 

Finally, the choice of a certain point within the hollow. This is a bit more 
difficult as computed quantities are far more sensitive to a variation within this 
line. As explained in Section ^, we take this as an indication of how far on can 
go with ones accuracy. We quote some representative numbers in Table ||, which 
are to be compare with /'(0) = -0.228 60 from the LPA. 

The computation of the eigenvalues is a quite easier problem, once the FP is 
determined (recall that now we are dealing with linear equations) . Here we shoot 
from the origin (with g'(0) fix) and fine-tune the initial condition h(0) and the 
eigenvalue A, by imposing that the eigenfunctions g((p) and h((p), together with 
their derivatives, remain bounded. (This somehow bold asymptotic conditions 
suffice.) Of course, one has to repeat the above analysis to be sure that the 
results are independent of the numerical method and, in particular, of the point 
where the asymptotic conditions are imposed. 
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